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Abstract 



In evolution equations for a complex amplitude, the phase obeys a much more intricate 
equation than the amplitude. Nevertheless, general methods should be applicable to both 
variables. On the example of the traveling wave reduction of the complex cubic-quintic 
Ginzburg-Landau equation (CGL5), we explain how to overcome the difficulties arising in 
two such methods: (i) the criterium that the sum of residues of an elliptic solution should be 
qq , zero, (ii) the construction of a first order differential equation admitting the given equation 

£N| ' as a differential consequence (subequation method). 

o ■ 

cn 

T^j-' ' 1 Introduction. Modulus vs. phase in amplitude equations 

The time evolution equation At + • • • = for a complex amplitude A(x, t) is usually, from 
physical requirements, invariant under an arbitrary shift of the phase <p = argA, in which M 
and ip denote the modulus and phase, A = Me lLp . As a consequence, in the coupled partial 
^ \ differential system for (M,(p), the variable (p only contributes by its derivatives. Then, under a 
reduction to an ordinary differential equation (ODE) such as the travelling wave reduction 



M -> M(£), p -> -tut + <p(£), £ = x- ct, (1) 

in the coupled ODE system for (M,ip = <p'), the highest derivation order for ip will be one less 
than the derivation order for M. Consequently, the ODE for ip obtained by the elimination of 
M will be much more complicated (by both its volume and its structure of singularities) than 
the ODE for M obtained by the elimination of ip. 

Let us take as an example the one-dimensional cubic-quintic complex Ginzburg-Landau equation 
(CGL5), 

iA t + P A xx + q\A\ 2 A + r\A\ A A-ijA = 0, (A,p,q,r) G C, pr / 0, 9?(r/p) ^ 0. (2) 



*To appear, Theoretical and Mathematical Physics. Solitions in 1+1 and 2+ldimensions. DS, KP and all that, 
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which depends on seven real parameters since 7 can be chosen real. For a summary of results 
on CGL5, see the reviews [HQS]. Its travelling wave reduction 



A(x,t) = ^/M(Oe 

/2 



H-ut + y>(0) 



, £ = x — ct, (c, oj, M, ip) G 72.. 

r ^2 , w-i 7 



— — - — =■ + »V - + V— - i — + -to' + -M + -M + '- = 0, 

2M 4M 2 ^ ^ M 2p M p p p p 



(3) 
(4) 



introduces two additional real constants (c, u), but the total number of real parameters (nine) 
can be lowered to seven by a translation of <p' and by noticing [TB] that the velocity c and 
the imaginary part of 1/p only contribute by their product. Indeed, denoting the eight real 
parameters as &r , 6^ , dy , ^ , , 5 % . Oy. , , 



. _ r . _ g . 1 

p p p 



i + — L {2s i + is r ), 5 

p 4 



and performing the translation 



/ cs r 1 



the system only depends on the seven real parameters e r ,ei,d r ,di,g r ,gi,csi = k\. 
The coupled two-component system in the real variables (M, ip), 



M" M 



M' 9 o 
«i — - V + e r M 2 + d r M + 5i = 0, 



2M 4M 2 
M 

^ + ~ ^ + €iM2 + d%M ~ 9r = °' 
contains as highest derivatives M" and The following parity invariances of ([7]) 

g = 0: (M,^f)->(-Af,V,0, 

K i = 0: (M,^,£)^(M,-^-£), 

will be used later on. 



(6) 



(7) 



(8) 
(9) 



The elimination of tjj yields a one-line real third order second degree ODE for M [12 



2mG - G' 



G 



2M 2 (eiM 2 + <UM - g r y M 2 ' 

(G" - 2 Ki G) 2 - 4GM 2 ( ei M 2 + djM - c/ r ) 2 = 0, 

G = -MM" - ^M' 2 - ^-MM' + e r M 4 + d r M 3 + 9i M 2 , 



(10) 

(11) 
(12) 



while the elimination of M yields a third order fourth degree ODE for ip which contains 11053 
terms, whose dominant ones (in the sense of singularities, as developed in section [2|) are, 

efV 4 { [25e?V>V - 10e 2 (3e^' + 4e r V 2 )V^2 - 6^' 3 ef - 24e r e 2 V>'V 



+ ei (80e 2 - 112e 2 )^V - 32e r (5e 2 + e 2 )^ 6 



3e 2 ^' 2 + 12e r e^V 2 + 4(5e 2 + 2e z r )tp 



1.1.2 



„2n <4 



We 2 ipip2 + e 2 ^' 2 + lGej.ej'i/'V 2 + 4(5e 2 + e 2 )^ 4 | + subdominant terms = 0. (13) 
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As a by-product of the elimination process, the rational expression (|10p of ip in terms of 
M,M',M" is quite short, while the rational expression of M in terms of ip also involves the 
third derivative ip'" and is quite lengthy. This is why the phase of the complex amplitude A is 
qualified as a "slave" variable [16], because it allows one to easily compute (p from M but not 
vice versa. 

The purpose of this work is to explain on the above example how to overcome the difficulties 
created by the consideration of ip in two specific methods. 

The paper is organized as follows. In section El as a prerequisite study, we investigate the detailed 
structure of movable singularities of M and ip, some features of which had been previously 
overlooked. 

In section [3l we apply a first method successively to M and ip in order to build necessary 
conditions for M or ip to be elliptic. This proves easy for M and sets up additional questions 
for ip, whose solution is provided. 

In section we indicate how to correctly apply a second method (the subequation method 
|14| 15} 16]) to M and ip, in order to build a first order ODE sharing elliptic or degenerate elliptic 
solutions with the above third order ODEs (fTTj) and (fT3j) , 

In section [5j we simply present the solution of CGL5 in which the square modulus M is elliptic, 
whose obtention by the subequation method is detailed elsewhere [8]. 

2 Movable singularities of CGL5 

In this section, we enumerate all the movable poles of either M or ip, excluding those which 
represent a singular solution^ of either the third order ODE (|lip for M or the third order ODE 
(JSP for ip. 

The structure of singularities of ([2]) has been studied in [13 1 . A first type of singularity x = 
£ — £o - > is obtained by balancing the terms A xx and |j4| 4 ^4 in ([2]), 

A ~ A oX { - 1/2+ia \ A ~ A oX ( - 1/2 ' ia \ (-1/2 + ia)(-3/2 + ia)p + A*r = 0, (14) 
these algebraic branch points of A define four values of Aq and two of a, 

{eiAtf - 4e r 4 - 3 = 0, a = |< (15) 

A 2 =e 2] j 2er+ e ^ A , «= 2Cr 2 + £lA , A = V /4^T^, e\ = el = l. (16) 

At these singularities, the square modulus M = \A\ 2 displays four simple poles 

M ~ m oX ~\ m = Al, (17) 

and a direct investigation of the third order ODE (|11|) shows that M admits no other movable 
polefl 

1 A singular solution [3] of an ODE is any solution which cannot be obtained from the general solution. Such 
a solution must cancel an odd multiplicity factor of the discriminant of the ODE. For instance the ODE 
whose discriminant is GM 2 (aM 2 + diM — g r ) 2 , admits as singular solutions all those of G = 0, which must be 
rejected since they are not solutions of the system (J7J)- 

2 The Laurent series M — ±(3/(4e, )) 1,/2 x~ 1 + • • ■ must be discarded since it cancels G and thus represents a 
singular solution. 



3 



Let us now count the poles of ip by considering the real system ([7|). A first set of poles of ip 
arises from the simple poles of M, 



M = moXi 1 

■2 



/ "4 

1+ T + 



Xi + Oixi] 



2 

ejm _i eim Q 
ip = —^-Xi + — 5— «i + m 



4(l + e?mg) , 

4dj + 5eid r mQ — efdiiriQ 



4(1 + efmg 



+ o(xi), 



(18) 
(19) 



in which both invariances dHJ) — dSJ) require changing rriQ to —rriQ. This first set defines four different 
simple poles of ip when q is nonzero and only two when it vanishes (in which case M 2 obeys an 
algebraic equation admitting two double poles). 

A second set of poles of ip, not considered in [13J, arises from the movable simple zeroes of M, 
and this set is best computed from the system ([7]). This is 



1 



1 



M M 

i 3 -l 
^=2 X2 

1 



X2 1 



1 + M lX2 + \mI + K[Mi - J -g r + | ft j xl + 0{xl) 

l + (m + Mi) X 2 + \m\ + Ik-Mi + \gi - ysr + \k\ \ x \ 



(20) 



o ^ (5r +J5i)Ki + 



3i«f (3di-jd r )M (lljK 2 +Ag r + Ajg i )M 1 



+ 



+ 3j Ki M 1 2 +jM 1 3 } X | + 0(x 



(21) 



in which Mo and Mi are arbitrary constants, and j is any square root of —1. Again, both 
invariances (|8])-(j9|) require changing Mq to —Mo, with the additional requirement Mi = when 
Ki = 0. This second set defines 2A^ simple zeroes of M, and either 2A^ (when q ^ 0) or (when 
q = 0) simple poles of ip, with N an undetermined integer. 

Remark 1 In the nonlinear Schrddinger equation (p real, q real, r = 7 = 0), the first integral 
Mifi = constant implies that the only poles of ip arise from the zeroes of M. 



We have searched for possible additional poles of ip by directly investigating the third order 
ODE ()13p for ip{Cj- One thus finds three kinds of families of movable simple poles, namely the 
two previous kinds (|19p and (|2ip plus the following third kind comprising two families 

, 4(e 2 + 5e 2 )pg - 16e r eip + 21e 2 = 0. (22) 

However, this series cancels a factor of the discriminant of (|13|) having an odd multiplicity, 
therefore it is a singular solution of ()13p and it must be discarded because it is not a solution of 
the system ([7|). 

Therefore another distinction between M and ip is the difference of complexity of their singularity 
structure: exactly 4 Laurent series for M, an unknown number of series for ip. 



i> = pox~ 



l + j + 0(x) 



Remark 2 Out of the two special values q = and = allowing an invariance in the 
differential system, see only the value q = is involved in the structure of singularities. 

The value k\ = will show up in next sections^ and\5i 
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3 On the characterization of nondegenerate elliptic solutions 



This is a classical result that, inside a fundamental domain, the sum of the residues of any 
elliptic function at its poles is equal to zero. This allowed Hone [11 J to take advantage of the 
Laurent series to generate the following necessary conditions for the solution u of an ODE to be 
nondegenerate elliptic, 

Vj£N,Vk£N: CJ k = residue ((u^)^ =0, (23) 

Laurent series 

in which the sum extends to any subset of the set of Laurent series, and thus to isolate those 
parameters for which the ODE might have a nondegenerate elliptic solution. 

If one assumes that M is elliptic, it follows from (jlOp that ip is also elliptic, so one can a priori 
use the Laurent series of either M or i/j to generate necessary conditions for M and ip to be 
elliptic. However, as already noticed in section [51 the situation is much more complicated with tp 
than with M. Indeed, some Laurent series of ip (those near ^2) depend on arbitrary constants, 
therefore one must first compute the number of different series (|21[) . then one must also solve 
the generated residues conditions for these extra arbitrary constants. Let us apply the criterium 
successively to M and to ip. 

3.1 Criterium of residues applied to M 

The variable M presents two advantages over t/i: it admits exactly four Laurent series (|18|) . 
and no arbitrary coefficient enters these series. The sum (|23p can include one, two, three or 

four Laurent series. However, with one or three series, the condition (|23|) applied to (m(°>) 
generates m,Q = 0, which is forbidden. With two series, because of the invariance mo — > —tuq 
of (|15p . the generated conditions are identical to those with four series, therefore only one case 
remains to study, that with four series in the sum (|23p . After computing the first seven terms 
of each series ([15]). one generates ten conditions Cj^ = 0, jk = 01, 02, 03, 04, 05, 06, 07, 12, 13, 
22, 

C5Y = 0, C$f = K ier = 0, Cof = Ki (l9e r di(ei 2 + 16e£) - eid r (9ei 2 + 6e*)) = 0, 

Cq4 = K,[P(e r ,ei,d r ,di, g r , gi, k?) = 0, = ^iP(e r , e«, d r , di, g r , gi, nf) = 0, ••• (24) 

in which the two P's are polynomials containing respectively 18 and 16 terms. In the case k- x 7^ 0, 
the generated constraints are 

«j / : Cf£ = e r = 0, C^i = d r = 0, = 16 gi + 3k 2 = 0, Cf/ = d l = 0. (25) 

We will see in section O that these necessary conditions are sufficient because there does exist 
an elliptic solution when they are obeyed. 

In the case k\ = 0, the parity invariance ([9|) selects very few nonidentically zero expressions 
Cj'l = 0, the first few ones being: C^f = P2, 2,5, 5,15,15 (dr , 9i, d r , di, e r , e^jl (140 terms, requires 6 
terms in the series), C3I (not computed, 8 terms in the series), C-ff (not computed, 10 terms 
in the series), etc, thus requiring the computation of many terms in the series to generate the 
necessary conditions. 

Fortunately, the unpleasant feature of having to deal with the second set of poles X2, whose 
number is unknown and which adds arbitrary coefficients to the system of necessary conditions, 

3 The indices denote the degree of P in its arguments. 
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can be avoided by associating to M a "subdominant" contribution of the "slave variable" ifi 
in the following way. Any product which is holomorphic near X2 enjoys 

the same properties as M, i.e. : to have only four Laurent series at xi, to introduce no extra 
arbitrary coefficients. 

The case = can then be settled easily, the simplest necessary conditions being, 

' i^M : 3ei(2e 2 + 3e 2 )d r - e r (8e 2 + lle 2 )di = 0, 

e r (8e 2 + lle 2 ) 

solved as d r — 



3ei(2 e 2 + 3e 



2\ ' 



ipM 2 : 32(2e; + 3ef) e^e^ + 4e r # r ) - (16e; + 19ef)(4e; + 3ef)e r df = 0, 

. , 4e r (16e 2 + 19ef)(4e 2 + 3e?) j2 

solved as % = q r H ^ — £ ott^ e r a„- , 

yj 3e* 96e 2 (2e 2 + 3e 2 ) 2 

V>M 5 : dje r P 2 ,4,i2,i4(9r,di,e r ,ei) = 0, 

V> 3 M 3 : £^2,4,14,16 (SV, d», e r , ej) = 0, 

k ipM 6 : e r P 3) 6,i8,2l(5r,^,e r ,e i ) = 0. 

With only seven terms in each series, this defines three sets of necessary conditions 

led 2 



Ki 
Ki 



0, e r 

0, Ck 
0, di 



3s 29e _ d 2 

"TTCj) Q>r — ~rz~di, 9r — ~Z ; 9i 
2 15 osi 



d r 

d r 



9i 
9i 



9r 



0. 
0. 



(26) 



the last set being the particular case Ki = of the conditions 
could well be refined by using more terms in each series. 



1 = 1, (27) 

(28) 
(29) 

The first two sets l|57 |> - (|25p 



3.2 Criterium of residues applied to ijj 

The motivation of this study comes from a previous application of the criterium to ip, in which 
Vernov [18] found the set of necessary conditions 



dip — dj 



9i9r 



0. 



(30) 



which only detects the subcase Ki = of the elliptic solution (@7|). This result was achieved 
from the Laurent series of tp 1 , j = 1,2,3,4, but only from the series (|19p near xii discarding 
the second set of poles (f2~T|) near \2- Let us show in this subsection that, if one takes account 

of both sets of poles of ip, the single consideration of the variables 
the correct results. 



allows one to recover 



Since an elliptic function possesses as many zeroes as poles, M must have exactly four simple 
zeroes ()20p . The corresponding set of poles of ip, necessarily simple, are then the following: 



1. - (q 0) 4 poles xi pl us 4 poles X2] 

2. - (q = 0) 2 poles xi pl us 2 poles X2- 



Let us denote the triplets (Mq, M\, j) in (j2~T|) as either the four sets (-Mo,i±, M\ t x±, ±i), 
(M],2±) Ml,2±j (case g ^ 0), or the two sets (Mq } ±, Mx t ±, ±i)(case q = 0). 

Let us denote Cj^f the sum of the four residues associated to the four poles xi 021), and 
the sum of two residues associated to two of the poles X2 ([21]) with opposite values of j. 
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The first set of sum of residues evaluates to |18] 



°oi 



4e r 



e; 



C 



02 



a 



,ipp 



03 



6— ^ + 4-%g r + 

ef 



9e| + 32e 2 2 520e 2 e^i + 256e^ + 303ef e 2 + 9ef j2 



64e 3 (e 2 + e 2 ) 2 

132ef + 229e 2 e 2 + 112 e 4) Q e 2 (64e^ + 135e 2 e 2 + 81e|) ^ 2 

C/-?' di<Y> I O « j / n O \ o (X^> 5 



12 



9i 



IQeHel + e 2 ) 2 ^ T ^ ° 64e,(e 2 + e 2 ) 2 
3 o H2e| + 23e 2 e 2 - 224e^S j2 



e r ° 2 
a 16 



c 2^2 



+ 



96e 2 (e 2 + e_ 

132 e 4 + 229e 2 e 2 + 112e^ 19ef - 123e 2 e 2 - 232ef 

Cr> CL<i (jb<r> | 



32(e 2 + e 



2\2 



32 ei (e 2 + ef 



2\2 



and the second set to (denoting M ls = M 1)+ + Ml _, M 1(i = M 1)+ - Mi 



^01 



u > o 02 



Mi 



<i>z 



' °03 



-5r 



3«M 



Id 



(Mi s + 2Ki), 



,1/jZ 
12 



°04 



2 ' Ud 4 
13 41 9 3 

-Mf s + —K[Mf s + - Kl 2 M ls + -«j 3 + -(2«i + M ls )M 2 d 

- ig r M ld + 5i (2 Ki + M is ) + -j(M ,+ + M ,_) + ^(Mo,_ 

A Ml 3 s + ^ Ki M 2 s + ^k 2 M 1s + + ^(2/q + Mi s )M 2 d 

5 
" 4 



Mn 



16 



16 



5 ig r M ld + + M ls ) + ^(M ,+ + M ,_) + ^(M 0i+ - M _), 

4 16 16 



(31) 



(32) 



The sets of conditions to be solved are then (the arguments of C - 
unknowns) 



fZ, 
k \ 



describe the additional 



q + 0, N = 2 : C ff + Cj* (M ,i± , M M± , ±i) + C]£ (M 0)2± , M 1)2± , ±») = 0, (33) 



ipZ, 



0, iV 



1: i<^r + C^(^o, ±! M 1)± ,±i) = 0. 



ylpZ 



(34) 



The first set of conditions ()33p contains too many unknowns to be solved when only seven terms 
in each series are considered, we leave the computation to the interested reader. 

For the second set of conditions (f34"|) . since q is zero, the Laurent series (f20j) - (f2"Tj) must possess 
the invariance ([5]), 

q = 0: (M,^,X,M ,M 1 ,j)^(-M,^,X,-M ,M 1 ,j), (35) 
and, when k; is also zero, the additional invariance 

q = 0, Ki = Q: Mi = 0, (M,il>,x,M Q ,j)->(M,-i/>,-x,-Mo,j). (36) 
This second set of conditions evaluates to 

ip 1 : e r = 0, 
V> 2 : 2M ls + Ki = 0, 
V> 3 : 8g r + 9 Ki M ld = 0, 
< V 5 : 128e i (M 2 + + M 2 _)-9^ i M ld (43K 2 + 272 5i ) = 0, (37) 
V> 4 : Ki(3K 2 + 16 5i ) = 0, 

^' 3 M 3 : M ld (128iei(M 2 + - M 2 _) + 1296c/iK 2 + 279/sf) = 0, 
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and only admits the two solutions 

3 2 ,,2 «i 2 (16op±9i«?) w (-9k, 2 ± Wig r ) . , 

«i ^ : e r = 0, 5i = --«f , M 2 ± = -^-f- ^, M 1± = ^ l — ^1, 38 

lb b4ej 3o/-q 

Ki = : e r = 0, 5 r = 0, M 1± = 0, M 2 + + M 2 _ = 0, (39) 
therefore the subcase K; 7^ does succeed to generate the desired constraints (|25p . 

Remark 3 The denominator Ki in A38\) corresponds to the factor k; of the M' 4 term in \4&\) 
and, when Ki vanishes, the singularities \2 do not exist any more in UM , see [JW . 



When one now requires both sets of necessary conditions to hold true (the set (|25p. (|27|) — f)29p 
generated in subsection 13 .1\ and the above set ([38]) - ([39]) ), then only two possibilities remain for 
an elliptic solution M to exist: either (|25p . which does define a nondegenerate elliptic solution, 
or e r = d r = d{ = g r = g% = 0, which defines the four rational solutions M = mo/(£ — £o),ip = 
( ei m 2 /2)/(e-Co). 



4 Method to determine all elliptic and degenerate elliptic solu- 
tions 

Consider an JV-th order autonomous algebraic ODE (|40p 

E( U ( N \...,u',u) = 0, ' = i-, (40) 

ax 

admitting at least one Laurent series 

+00 

u = x p ^u j x' , X = x - x . (41) 
3=0 

There exists an algorithm [14J to find in closed form all its elliptic or degenerate elliptic solutions. 
Its successive steps are [SI [5] : 

1. Find the analytic structure of movable singularities (e.g., 4 families of simple poles, 2 of 
double poles). For each subset of families (e.g. 2 families of simple poles) deduce the 
elliptic orders m, n (e.g. m = 2, n = 4) of u, u' and perform the next steps. 

2. Compute slightly more than (m + l) 2 terms in the Laurent series. 

3. Define the first order m-th degree subequation F(u,vf) = (it contains at most (m + l) 2 
coefficients aj^), 

m 2m— 2k 

F(u,u) = ^2 a 3\kU j u =0, a , m /0. (42) 

k=o j=o 

According to classical results of Briot-Bouquet and Painleve (see details in [6]), any el- 
liptic or degenerate elliptic solution of (Ii0|) must obey such an ODE, which is called 
"subequation" of (|3D|) because it is required (in next step) to admit (jlOj) as a differential 
consequence. 
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4. Require each Laurent series (jiTj) to obey F(u, u') = 0, 

F = X m ^ \J2 FjX j + 0(x J+1 ) j , Vj : F, = 0. (43) 

and solve this linear overdetermined system for aj^. 

5. Integrate each resulting ODE F(u, u') = 0. 

The structure of singularities of M and ip has been established in section [21 and the result of 
concern to us is: the movable poles of M and ip are all simple, and the number of distinct 
Laurent series at these simple poles is equal to: 4 in the case of M, 4 + 2N in the case of ip or 
M' /M if q 7^ 0, and 2 + N in the case of ip or M'/M if q = 0, with iV an undetermined integer. 

Let us start with M. In order to find all elliptic and degenerate elliptic solutions M by the 
subequation method, at step 1 the various subsets of families to be considered are made of 
one, two, three or four series (|18p . defining subequations F = whose degrees in (M',M) 
are respectively (1,2), (2,4), (3,6), (4,8). The computation presents no other difficulties than 
technical ones and its detailed results can be found in [8]. The main new result is a nondegenerate 
elliptic solution presented in section [5j 

Applying the subequation method with ip presents several difficulties. 

1. The main one is not to forget the movable singularities of type \2 (movable simple zeroes 
(|20p for M, movable simple poles (|2ip for ip). This is the reason why a previous inves- 
tigation [18| could only find a particular case of the elliptic solution presented in section 

m 

2. The number of Laurent series (|2ip for ip is undetermined, thus failing to set an upper 
bound to the degree of the subequation F. This point has already been settled in section 
13.21 where we have established that the number of distinct Laurent series ip is either eight 
(q / 0, four series near xi plus four series near \2) or four (q = 0, two series near xi plus 
two series near X2)- 

3. The arbitrary coefficients in the series (|2ip must be determined by the subequation method 
and therefore require the computation of many more terms in each Laurent series. 

4. Another difficulty, undetectable a priori, is the infinite value of M\± at k\ = in (|38p . 
Setting K{ = (for some reason) while looking for subequations being obeyed by four series 
will fail. In such a case (q = and Ki = 0), ip admits two distinct series near xi an d n ° 
series near X2, and the correct assumption for F is 

F = ip' 2 + (a i + ax\ip + a 2 iip 2 )ip' + a o + aio^ + a20^ 2 + a30^ 3 + a 4 oip 4 = 0, (44) 

which is further reduced by a classical theorem (no first degree term |2J §181] in ip' since 
ip is assumed elliptic) and by the invariance ([9]) to the canonical Briot-Bouquet type 

F = ip' 2 + a 00 + a 20 ip 2 + a m ip A = 0, (45) 

for which a solution was indeed found [IS], see (|4"9"j) . 

For all those reasons, the choice of M is by far the best one for amplitude equations such as 
CGL3, CGL5 or others. 
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5 The elliptic solution of CGL5 



By "elliptic solution of CGL5" , we mean a solution of CGL5 in which the square modulus M of 
the traveling wave reduction is elliptic. 

The subequation method yields a unique genus-one subequation for M [8|, and it requires exactly 
the constraints (|25|) . Introducing the shorthand notation, 



61 = 48' 60 = 36' (46) 



this fourth degree subequation is 



M' 4 - 2 Ki MM' 3 + — ei M' 2 (eiM 2 - 12e ) + 

ei ef 

=,2 



+ (288e 2 + 24 ei e M 2 - e?M 4 ) - M 2 (e;M 2 - 48e ) 3 = 



(47) 



Because of the previously mentioned numerous difficulties with tp, the corresponding subequation 
for iff is not determined by the subequation method, but by elimination with the correspondence 
(fTOj) . resulting in 

' /qV/ 4 - 4/q?// 3 + 24e o) 

+ Sip' 2 (-Ki(27e\ - 324eg) + 1440eie V + 27K i e 1 if) 2 + 16e ^ 3 + g«i^ 4 J 
1 32 

+ 16( - -/qV 8 - ye ^ 7 - 26Kiei^ 6 - 1632eie ^ 5 - (477e 2 + 552eg) V 4 ( 48 ) 
- 288 (l65e 2 + 4e§) V 3 + «i (2106e 2 - 31320eg) ^ 2 
+ 2 7 3 6 (e\ - 4eg) e^Qip + 243 (-9ef + 56e 2 e§ - 144e<j) ) = 0. 



The singularity k x = already uncovered in (j38|) is displayed as the factor k\ in front of the ip /4: 
term in (ji"8]) . with the consequence that the degree of subequation (|4"8"j) drops from four to two 
when K[ = 0. One then recovers the result of Vernov |18| . 



0, e r = 0, 9l = 0, „ = (,,< ~f 3e,M 2 - 4*) = 0, 

9V' - 12^ 4 - 9r = 0, 



(49) 



in which both subequations for M and ip belong to the list of five canonical equations of Briot 
and Bouquet. 

Full details on the integration of (|47p and (|48|) can be found in [5]. The final result for the 
complex amplitude A is, 



— Ibjt + I ^ 

y Ki : A = constant e 2 P H(£, , 0)(- 1+ ^/ 2 H(£, O)^ 1 -*^/ 2 , 



(50) 



in which H(£, q, k) is the element simple defined by Hermite [101 vol. II, p. 506] for integrating 
the Lame equation, 



(51) 
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and the fixed constants £± in (|50|) are defined by 



P(C±, G 2 , G 3 ) = -2ei ± iV3(3ei + 4ie ), p'(£| , 
G 2 = 12(13e? + 16eg), G 3 = 8(35e 2 + 48eg)ei, 
p' 2 = 4(p + 2ei)(p 2 - 2e x p - 35e 2 - 48eg). 



G2, G3) 



3jFn/3 
2 



(3ei + 4ieo), 



(52) 



Numerical simulations with periodic boundary conditions |15(, Fig. 4] do display solutions M hav- 
ing a real period (similar features are observed in CGL3 O Fig. 7]), these could well correspond 
to the present elliptic solution. 

6 Conclusion 

The traps described in this article should be kept in mind when looking for all the elliptic or 
degenerate elliptic solutions of other amplitude equations, such as the complex Swift-Hohenberg 
equation [17 1 . 
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